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ABSTRACT 

We perform a frequentist analysis of the bispectrum of WMAP first year data. We 
find clear signal domination up to t sa 200, with overall consistency with Gaussianity, 
except for the following features. There is a flat patch (i.e. a low \ 2 region) in the 
same-£ components of the bispectrum spanning the range I = 32 — 62; this may he 
interpreted as ruling out Gaussianity at the 99.6% confidence level. There is also an 
asymmetry between the North and South inter-/! bispectrum components at the 99% 
confidence level. The preferred asymmetry axis correlates well with the (£, b) = (57, 10) 
direction quoted in the literature for asymmetries in the power spectrum and three- 
point correlation function. However our analysis of the quadrupolc (its bispectrum 
and principal axes) fails to make contact with previously claimed anomalies. 

Key words: Cosmic microwave background - Gaussianity tests. 



1 INTRODUCTION 

The hypothesis of statistical Gaussianity plays a central role 
in the standard theory of structure formation. Unsurpris- 
ingly, a large amount of effort has gone into checking this 
assumption, with the cosmic microwave background (CMB) 
anisotropy being its cleanest probe. The early work on 
COBE-DMR data was punctuated by the first findings that 
we might not live in the best of all possible worlds. The data 
passe d standard, but n ot very discriminative Gaussianity 
tests jKogut et a]|[l9 96). Upo n closer inspection, however, 
anomalies were discovered <|F erreira 1 Maguei io fc Gorskil 
Il998l: iNovikov. Feldman and Shandarinl ll99Sl) . For exam- 
ple, the bispectrum analysis o f the 4-year dataset revealed 
evidence against Gau ssianity jFerreira. Magueiio fc Gorskil 
ll99StlMagueiidl200Ch to a high significance level. 

The possibility remained that these non- Gaussianity 
"detections" might be due to undocumented systematic 
errors. This hope was partly satisfied by the discov- 
ery of subtle errors in the map-making pr ocedure of 
COBE-DMR ijBandav. Zaroubi fc Gorskil Il999h . Pixeliza- 
tion artifacts and further systematic errors in COBE 
were later fo und, with the aid of the higher quality 
WMAP data JMagueiio fc Medeirosl |2003l) . Indeed all the 
early work was plagued by high noise levels, poor resolu- 
tion an d the threat of mul tiple systematics. The WMAP 
proiect iBennett fc all l2003aTl has opened up the doors for 
the first convincing evaluation of the hypothesis of Gaus- 
sianity. 
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worry- 

ingly, it appears that the data displays a North-South 
asymmetry in the power spectrum and three-point correla- 
tion function. This questions both the statistical isotropy 
and Gaussianity of the data. The spread in power spectrum 
and 3-point function is sensitive to the fourth and sixth 
momenta of the distribution. 

In this paper we perform a frequentist analysis of the 
W MAP bispectrum. Th is is to be contrasted with the work 
of iKomatsu fc all <2003fl . where the WMAP bispectrum is 
analysed with reference to a particular non-Gaussian model, 
parameterised by the non-linear coupling parameter /jvl. 
Such work neglects bispectrum non-Gaussianity that might 
arise in any other theory; or indeed that might be due to spu- 
rious systematics. Our exp loratory approach is t o be seen as 
complementary to that bv lKomatsu fc al (|2,003). In partic- 
ular we shall be concerned with correlating our bispectrum 
studies and the WMAP anomalies found by other groups, 
namely claims for North-South irregularities. 

Our paper is organised as follows. In Section |5] we es- 
tablish our notation and bispectrum definitions. We intro- 
duce same-£ and inter-/! normalised bispectrum measures, 
denoted by It and Ji respectively. These are sensitive to 
phase correlations within a single multipole and between 
adjacent multipoles. Their normalisation ensures that for a 
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Gaussian process they are statist ically independent from the 
power spectrum jMagueiidll995ft . 

In SectionEJwe then evaluate Ie and Ji for WMAP and 
compare results with Monte-Carlo simulations in which a 
Gaussian signal is subject to the same noise and beam char- 
acteristics as the real data. We also discuss various galac- 
tic masks. We find a low \ 2 region in the It in the range 
£ — 32 — 62. Depending on the prior, this may be inter- 
preted as ruling out Gaussianity at the 99.6% confidence 
level (we relegate to Section HJ a discussion of possible selec- 
tion effects). 

In Section 2] we re-examine claims for North-South 
asymmetries in CMB from the point of view of the bis- 
pectrum. We find asymmetries in the inter-^ components 
of the bispectrum which are not only statistically significant 
(99% confidence level, even when all possible selection ef- 
fects are taken into account), but also correlates well with 
the (£, b) = (57, 10) direction quoted in the literature. Our 
detection is independent from, but reinforces si milar anoma- 
lies fou nd for the 3-point correlation function jEriksen fc all 
l2004aft . 

However, our analysis of the quadrupole fails to make 
contact with previously claimed weird features. In Section|S] 
we provide a geometrical interpretation of the quadrupole 
where the bispectrum may be seen as a shape factor. In 
addition we identify three orthogonal principal directions. 
None of these show anything abnormal or correlate with 
suspicious directions. 

In Section we conclude with a summary of our find- 
ings. 



2 THE BISPECTRUM 

We start by reviewing some results and definitions pertain- 
ing to the bispectrum. Given a full-sky map, ^£(n), this 
may be expanded into Spherical Harmonic functions: 



— (n) = 2J STt = 2J a em Y em (n) 



(1) 



Each multipole component STi contains fluctuations with 
a characteristic angular scale inversely proportional to t. 
The coefficients atm are not rotationally invariant but may 
be combined into rotationally invarian t multilinear forms 
fsee lMagueiio. Ferreira fc Gorskil (1 1998ft for a possible algo- 
rithm). The only bilinear invariant is the commonly used 
angular power spectrum, given by C e = J2 m \ a c™\ 2 ■ 
This gives a measure of the overall intensity of each mul- 
tipole. The most general cubic invariant is the bispectrum, 
and is given by 
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where the (...) is the Wigner 3 J symbol. The proportion- 
ality constant is chosen so as to enforce a roughly constant 
cosmic variance. There are various alternative ways in which 
the bispectrum may be computed; for instance 



Bt 1 i 2 £ 3 oc / dflST^ 5Ti 2 8T1, 



(3) 



with s imilar constructions for higher orde r mo - 
ments dSpergel fc Goldberd Il999t iKomatsu fcall 12002ft . 

This may be computationally advantageous. 

Selection rules require that £\ + £2 + £3 be even, and 
we impose £i ^ 2. So for even £ the choice l\ = £2 = 
£3 = £ leads to the "single-l" bispectrum Be — Beee 
iFerreira. Maeueiio fc G6rskill99lift . Other bispectrum com- 
ponents are sensitive to correlations between different scales. 
The simplest chain of correl ators is therefore At = Be-ne+i 
- the "inter-i"' bispectrum dMagueiid l2000*l . Other compo- 
nents, involving more distan t multipoles, may be considered 
JSandvik fc MagueiidExIlft but they are very likely to be 
dominated by noise; it has been observed that non-Gaussian 
inter-scale correlations decay with £ separation. 

In this paper we shall consider ratios 



B, 



N3/2 



and 



T 3 - 
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(4) 



(3) 



These quantities are invariant under rotations and parity; 
they are also dimensionless and therefore less dependent 
upon the power s pectrum. Indee d for a Gaussian process 
it may be proved <lMagueiidll995ft that it is this normalised 
- and not the un-normalized - bispectrum that is statisti- 
cally independent from the power spectrum. This matter will 
be very important later when we compare previously found 
North/South power spectrum asymmetries with asymme- 
tries found in the normalised bispectrum. 



3 THE WMAP BISPECTRUM 



The WMAP mission feennett fc all l2003aft produced full 
sky CMB maps from ten differencing assemblies (DAs); 
four in the W band (94 GHz), two V band (61 GHz), 
two Q band (41 GHz), one Ka band (33 GHz), and one 
K band (23 GHz). The K and Ka bands are dominated 
by galactic emission and therefore neglected for cosmo- 
logical analysis. We use the "foreground cleaned" maps 1 , 
where the Galactic foreground signal, consisting of syn- 
chrotron, free-free, and dust emission, was removed using 
th e 3-band, 5-parameter template fitting method described 
in lBennett fc al (2003b). All maps are rendered usi ng the 
HEALPix 2 package (iGorski. Hivon fc Wandeltlll99Sft . with 
pixelization level nside=512. 

We examined the bispectrum of an "inverse- noise- 
squared" coadded map based on the Q, V and W bands, 
that is a map: 



T(n) 



E^ 3 V-! 



(6) 



1 Publicly available from |http: / /lajnbda.gsfc.nasa.gov| 

2 http:/ /www. eso.org/science/healpix/index. html 
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where 

Npix 
n=l 



erf (n) = ag ti /Nobsi(n) 



where T; is the sky map for the DA i with the foreground 
galactic signal subtracted, and CTq i is * ne no ise per obser- 



vation for DA i, whose values are given by Bennett & a 
( 2003aj) . We apply the KpO mask to this map l Bennett fc a 
I2003VJ) to cut possible remaining galactic contamination 



near the disk, and then we remove any residual monopole. 

In Figs.0and|5|we plot the if and jf bispectrum com- 
ponents, for I = 2 — 450, and £ = 4 — 450 respectively, for the 
WMAP first year coadded map with a KpO mask. We com- 
pare these results with those from Monte Carlo simulations 
of Gaussian maps with the same noise and beam charac- 
teristics as the WMAP instrument, and subject to the KpO 
mask. From 5,000 sky realizations we inferred the 1, 2, and 
3-sigma variance lines, also plotted in Figs. 0and|21 

Before we apply a goodness of fit statistic to these re- 
sults it is important to select a signal dominated I range. 
This should be estimated from real data; any other algo- 
rithms risk including noise-dominated multipoles into the 
analysis, thus ending up testing the Gaussianity of noise. 
Such a direct estimate may be obtained by evaluating if 
and Jf for the different channels. The noise in each of these 
bands is uncorrelated, so the impact of the noise upon the 
statistics may be inferred from the dispersion from channel 
to channel. As can be seen in Figs. |3] and |3] the bispectrum 
is signal dominated up to £ w 200. 

Bearing this in mind we may now compute the reduced 
chi-squared for a selection of signal dominated sections in £ 
space. Table^records the fraction of simulations that found 
a lower \ 2 value for the indicated I range (note that the I 
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Figure 3. The WMAP single-^ bispectrum for the channels that 
most contribute to the coadded map (Ql, Q2, VI and V2). 



range can start at i — 4 for Ji and at £ = 2 for I(). We find 
the WMAP Jf consistent with Gaussianity. However, for the 
if we observe a long "flat patch" : a region where all points 
are within 1 sigma. Our "flat patch" covers points £ — 32—62 
inclusively, as can be seen in Fig.0 For this range we would 
expect a low x 2 value, and at 0.329 it is smaller than 99.6% 
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Figure 4. The WMAP inter-i? bispectrum for the channels that 
most contribute to the coadded map (Ql, Q2, VI and V2). 



lutes our "flat patch" . For both maps the "flat patch" does 
not diminish as we increase the masked region, implying that 
the "flat patch" is not due to galactic emission. 

We conclude that the the normalised inter-^ and single- 
£ bispectrum from the whole sky are consistent with Gaus- 
sianity. 



4 HEMISPHERE ASYMMETRIES IN THE 
BISPECTRUM 

Recently there have be en numerous detections of a s ymme - 
try in the CMB (e. g.. ICopi. Huterer fc Starkmanl (1200: 
IColes fc all j2003l); lEriksen fc all _ J2004bl) : lHansen fc a : 



Eriksen fc all 



Eriksen fc all 



11200411: lHansen. Banday fc Gorskil J2004I) : 
l|2004a[) : ISchwarz fc all i2004h ). ^Specifically. 
( 2004a|) found a significant lack of power and a surprisingly 
featureless pseudo collapsed 3-point correlation function, 
C^ 3 '(0), in the Northern ecliptic hemisphere. The asymme- 
try in the power spectrum was maximised for a North pole at 
approximately {£, b) = (57, 10) 3 . The non- normalised bispec- 
trum is related to the Legendre transform of C' 3 - 1 (0) . There- 
fore by observing the if and Jf on different hemispheres we 
can probe any £ dependence of the C' 3 ' (9) asymmetry. We 
emphasise that these will also be i ndependen t of th e power 
spectrum for a Gaussian field (see iMagueiiol (1 19951) for an 
anlytical proof for £ — 2) . 



Table 1. The fraction of our simulations that found a lower x 2 



1 range 


j3 


I 3 


2(4)-200 


0.75 


0.62 


2(4)-100 


0.48 


0.58 


2(4)-50 


0.466 


0.617 


30-70 


0.239 


0.033 


32-62 


0.141 


0.004 


50-100 


0.409 


0.437 


76-126 


0.758 


0.642 


100-150 


0.854 


0.445 


126-176 


0.79 


0.82 


150-200 


0.746 


0.617 



of those from simulations for this range. Critically we need 
to take account of selection effects, and so by letting each 
simulation pick out its longest "flat patch" within £ = 2—300 
we can account for the fact that we have honed into the most 
significant region. We find that a "flat patch" of this length 
or longer appears in 8.8% of the simulations. Further, a \ 2 
of 0.329 is not at all significant for a "flat patch" the same 
length as ours: 67.7% having an even lower x 2 ■ Of course 
the last conclusion depends on whether or not a prior exists 
which makes the range £ = 32 — 62 special. If such a prior 
exists then we have rejected Gaussianity at a high level of 
confidence. We shall return to this matter in Section |S| 

It is important to evaluate the impact of the galac- 
tic cut on these con clusions. We note that the ILC map 
jBennett fc al [2003a) also displays this "flat patch", but 
only when the central galactic area is masked. The ILC map 
does not completely avoid foreground emissions, and so the 
galactic region will display non-Gaussian signals. When the 
ILC map is not masked, it is this galactic emission that di- 



4.1 Inter-£ normalised bispectrum 

We calculate the Jf on the Northern and Southern hemi- 
spheres for a North pole at (£, b) — (57,10). In Fig. we 
plot them against each other, and against the whole sky re- 
sults. We also plot the 1 sigma variance lines for the whole 
sky (the spread is a little wider for a half-sky). We notice 
that the Northern results appear to take higher values than 
the Southern. For 34 of the 49 points the Northern result 
is higher, and when it is not the difference is small. We 
note that the individual \ 2 values do not reveal any sig- 
nificant departure from Gaussianity for either the Northern 
or Southern results. It is the difference between the hemi- 
spheres that is significant. 

We quantify this feature by measuring the sum of the 
differences: 

K = (7) 

t 

We calculate K for a North pole at (57,10), and at the eclip- 
tic North pole position (96,30), and at a further 53 positions 
uniformly distributed in the Northern galactic hemisphere. 
We infer the value of K for the opposite North pole posi- 
tions, in the Southern galactic hemisphere. We then perform 
Gaussian Monte Carlo simulations, as above, in order to 
evaluate the statistical significance of the observed WMAP 
asymmetry. For each simulation we allow the North pole to 
vary over the same positions and we record the maximum 
\K\ value seen in each simulation. 

3 Galactic coordinates: Longitude I S [-180(E), +180(W)] , Lat- 
itude b G [-90(5), +90(A% with {£, b) = (0, 0) towards the galac- 
tic centre. 
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Figure 5. The jf spectra: top panel, allsky (bold) and North- 
ern (dashed); middle panel, allsky (bold) and Southern (dashed); 
bottom panel, Northern (bold) and Southern (dashed). 
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Figure 6. Asymmetry in the j| as measured by \K\. At each 
North pole position we show the fraction of simulations with a 
lower l-Klmo! value, multiplied by the +/— sign of the Ky^MAP 
result. Results are also shown for North pole position (57,10) and 
its opposite (-123,-10), and for the ecliptic poles. 

Figure |5] plots, at the corresponding North pole posi- 
tion, the fraction of simulations with a lower \K\ max value 
than the \Kwmap\ value from that North pole. This fraction 
is then displayed with the +/— sign of this Kwmap result 
i.e. -1, and +1 indicates a North pole position returning 
higher asymmetry than all the simulations. We show the re- 
sults for the ranges I — 4—100 (larger circles), and £ = 4 — 50 
(smaller circles). 

In Table [5] we list the results for the three North pole 
positions that find the maximum value of K, for WMAP. 
We also show the fraction of simulations that return a lower 
| K | max value. 



Table 2. Measuring asymmetry in the Jf with K. The three 
North poles that maximise the asymmetry for WMAP, and the 
fraction of simulations with a lower value of \K\ max . 



1 = 4 - 150 


(57,10) 


0.808 


(63.9,20) 


0.707 


(85.1,20) 


0.689 


I = 4 - 100 


(57,10) 


0.994 


(63.9,20) 


0.946 


(40,0) 


0.832 


I = 4 - 50 


(63.9,20) 


0.988 


(57,10) 


0.922 


(40,0) 


0.916 



We observe significant asymmetry at the ~ 99% level 
for both £ ranges 4-100, and 4-50. The North pole positions 
that maximise this asymmetry differ slightly between the £ 
ranges, but b o th are c onsistent with the (57,10) position that 
lEriksen fc ail §004a) found maximises the angular power 
spectrum asymmetry. Furthermore, the pattern on the sky 
displayed in Figgis ve ry similar to that for the dis t ributi on 
of power, see Fig. 24 in (Hansen! Bandav fc Gorskil ll2004l) . 

Whether this anomaly is a sign of anisotropy, non- 
Gaussianity or both is partly rhetorical. An anisotropic 
ensemble may always be converted into an isotropic one 
by randomising the preferred axis. From the point of 
view of such an ensemble, what appears in a single 
sky as a sign of aniso t ropy is simply non-Gaussianity 
(see lFerreira fc Maeueiiol dl997l) for details). 

As previously mentioned, for a Gaussian random field 
the normalised bispectrum and the power spectrum are in- 
dependent quantities. In the analysis in this section we were 
particularly careful to avoid selection effects (we let Gaus- 
sian simulations choose their maximal asymmetry axis). 
Should one insist that the anomaly found is a fluke, how- 
ever, we should stress that it would be an independent 
fluke with respect to the a nomaly found in the power spec- 
trum feriksen fc alll2004al) . This makes it even more un- 
likely. We shall return to this matter in Section HJ 

4.2 Single-^ normalised bispectrum 

We calculate the single-^ normalised bispectrum, if, on the 
Northern and Southern hemispheres for a North pole at 
(£, b) =(57,10). In Fig. |7|we plot them against each other, 
and against the whole sky results. We also plot the 1 sigma 
error bars for the whole sky (the spread is very similar for 
a half-sky). The Northern results appear to fluctuate more 
widely than the Southern, and to investigate this further we 
calculate the \ 2 values for both hemispheres, for various £ 
ranges. 

In Table[!|]we show the fraction of simulations that find 
lower x 2 values, with our simulations also using a North pole 
at (57,10). The WMAP Northern \ 2 values are consistently 
high, while the Southern are low. For £ =2-100 the Southern 
X 2 is low at the 99% level. However, the significance of this 
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Figure 7. The if spectra: top panel, allsky (bold) and North- 
ern (dashed); middle panel, allsky (bold) and Southern (dashed); 
bottom panel, Northern (bold) and Southern (dashed). 



Table 3. Measuring asymmetry in the If with x 2 ■ The fraction 
of simulations which have a lower x 2 f° r the South, a lower x 2 
for the North, and a lower ratio of the two. 



1 range 


S 


N 


N/S 


2-200 


0.10 


0.88 


0.96 


2-150 


0.03 


0.91 


0.99 


2-100 


0.01 


0.88 


0.99 


32-62 


0.14 


0.88 


0.95 


2-50 


0.12 


0.95 


0.97 


2-40 


0.07 


0.83 


0.96 



asymmetry is not seen in these individual values from the 
North and the South, but in the ratio of the two: the fact 
that one is low while the other is high. Therefore, in Tabled 
we also record the fraction of simulations which find a lower 

ratio, MAX(^4,^f). 

x s x N 

Our if asymmetry appears significant at the 99% level 
for a North pole at (57,10). However, we must take account 
of our North pole selection as we may have pre-selected the 
North pole that maximises this asymmetry while we have 
not done the same for the simulations. As before, we vary 
the North pole position for the WMAP and for the sim- 
ulations through a further 53 positions that are uniformly 
distributed in the Northern galactic hemisphere (and infer 
the results for the positions on the Southern galactic hemi- 
sphere). At every North pole position we calculate Xjv/xIj 
and x% /XiVi an d f° r each simulation we record the maximum 
chi-squared ratio result from all these values. In Table 2] wc 
record the North pole positions which maximise the ratio 
for the WMAP data, and the fraction of simulations with 
lower maximum ratio values. We visualise this in Fig. |H| 



Table 4. Measuring asymmetry in the If with the North-South 
X 2 ratio. The fraction of simulations which, when varying the 
North pole, find a lower maximum ratio. 



1 = 2 - 150 


(57,10) 


0.581 


(340.5,20) 


0.180 


(149.0,20) 


0.108 


I = 2 - 100 


(180,0) 


0.683 


(57,10) 


0.527 


(149,20) 


0.431 


I = 2 - 50 


(0,20) 


0.940 


(340.5,20) 


0.587 


(21.3,20) 


0.557 
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Figure 8. Asymmetry in the if with the N-S x 2 ratio. Fraction 
of simulations with a lower maximum ratio. 



The significance of our asymmetry detection in the 
single-^ bispectrum has weakened greatly now that we have 
taken a selection effect into account. In fact the only hint of 
asymmetry that survives is at a level of 94%, for the range 
£ =2-50 with North pole at (70,0). This result is not signif- 
icant enough to rule out Gaussianity. In Fig. [S] we see that 
very few positions have a chi-squared ratio comparable to 
the simulations' maxima. Further, those that do, group to- 
gether towards the galactic centre. This is a very different 
pattern from that observed for the inter-^ bispectrum asym- 
metry, implying that any asymmetry here is unrelated. The 
grouping towards the galactic centre does raise some con- 
cern that the mask or foreground contamination may have 
a subtle effect on the single-^ bispectrum. 



4.3 3-point correlation function 

We briefly comment on the three-point corr elation func- 
tion, C ^(fl), and its asymmetry reported bv lEriksen fc all 
l|2004al) . We refer to this work for a direct definition (in 
position space) of C (3) (0). The C (3) (6>), is related to our 
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Figure 9. The threepoint correlation function, (9) in 

10 6 (fiK) s , from the Northern and Southern ecliptic hemispheres, 
with 1,2 and 3 sigma error bars from simulations. 



definition of the bispectrum Be 1 e 2 e 3 , by 
C {3> (9) = B ei e 2i3 Pi 3 (cos(9)) x 



(8) 



£i ii £3 \ (2^i + l) (26 + 1) (2^3 + 1) 




(47t) 



^2 Xw 2t3 Pi 3 {cos{0)) 



It is computationally advantageous to calculate C^\9) in 
Fourier space - through this summatio n - than directly over 
all points in the sky. Initially, we follow Eriksen fe~a 
and choose not to include the dipole and quadrupole terms, 
because of their well known anomalous behaviour; therefore 
we perform the summation in equation (|HJ for £ ^3. 

In Fig. |5] we plot the results from the Northern and 
Southern ecliptic hemispheres (scaled by their respective ar- 
eas), with 1,2 and 3 sigma error bars from simulations. As 
can be ap preciated from these p lots, we have confirmed the 
findings of lEriksen fcall (l2004al) : a featureless C (3) (6>) in the 
Northern ecliptic hemisphere. It is unlikely, however, that 
this result is simply related to the bispectrum asymmetry 
reported in Section f4.ll as we shall now argue. 

The three-point correlation function is a function in 
real space, while the bispectrum and power spectrum are in 
Fourier space. It is often not intuitively clear how a feature 
in one will manifest itself in the other. The \Xt t i 2 t 3 \ values 
decrease very quickly with increasing £, so it is the earliest 
values that contribute the most, and determine the over- 
all shape of C {3) (9). Including t ^8 terms returns C (3) (0) 
values approximately an order of magnitude less than when 
including I ^2 terms. Therefore, an effect that looks to prop- 
agate through a large range in real space, could be from a 



Table 5. Excluding lower I contributions in the calculation of 
C^(9). The fraction of simulations that found a lower \ 2 from 
the Northern and Southern ecliptic hemispheres. 



Min I 


X N 


4 


2 


0.066 


0.716 


3 


0.125 


0.912 


4 


0.068 


0.521 


5 


0.074 


0.442 


6 


0.000 


0.601 


7 


0.020 


0.724 


8 


0.009 


0.687 



small number of values at this low-£ end in Fourier space. 
Taking this into account could decrease the significance of 
such a detection of asymmetry. In light of current uncer- 
tainties about foreground con t amination in the low-l mul- 
tipole s (ISlosar fc SeliaklEx) J: [Bielewicz. Gorski fc Bandavl 



12004 lHansen. Bandav fc O6rskill2004h we investigate if the 
asymmetry survives when we exclude low-£ contributions. 

We repeat the calculation of C^{9) on the Northern 
and Southern hemispheres, and each time we increase the 
minimum I value allowed to contribute, i.e., start the sum- 
mation in equation JSJ from an increasing minimum i value. 
We do the same for our simulations. For each computation, 
we perform a reduced chi-squared analysis over the 9 range 
[0,70], and in Tabled we record the fraction of ~500 sim- 
ulations that find a lower \ 2 value than the Northern and 
Southern WMAP results. Our simulations also use the eclip- 
tic poles. We stress that we are not interested in the level 
of significance of this anomaly (we have not accounted for 
selection effects), but how the significance changes as we 
exclude low-£ terms. We find the asymmetry does not di- 
minish as we exclude the \ower-£ contributions, implying an 
asymmetry that is a widespread feature in £ space. 

We emphasise that this asymmetry is different to that 
explored above with if and jf , as here we are not indepen- 
dent of the power spectrum. 



5 THE QUADRUPOLE 

The geometrical meaning of the bispectrum becomes clearer 
by focusing on the quadrupole. Also, the low value of the 
WMAP quadrupole has inspired much speculation; our re- 
marks shed some light on this matter. 

In general a given multipole may be decomposed into 
the angular power spectrum Ct (the radius) and an axis 
(a direction), in 2£ + 1 dimensions. Alternatively the lat- 
ter may be represented as £ independent un i t vectors in 3 
dime nsions llCopi. Huterer fc Starkman||2003t ISchwarz fc all 
2004). These contain the remaining 2£ degrees of freedom 
(d.o.f); but given that there are 3 rotational d.o.f. in the 
choice of the coordinate frame (and that each multipole is 
an irreducible representation of the rotation group), such 
vectors cannot be rotationally invariant, i.e. they depend on 
the way the system of coordinates was set up. 

For instance the quadrupole has a total of 5 d.o.f. One is 
C2, so there is only one ot her invariant. The 2 3D axes result- 
ing from the procedure in lCopi. Huterer fc Starkmanl j2003ll 
contain all the remaining 4 degrees of freedom, so they mix 
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Table 6. The Quadrupole bispectrum and eigenvec- 
tors; and some "unusual directions". 



Dataset 








V 2 




V3 








I 


b 


I 


6 


I 


b 


Coadd, KpO 


-.458 


48.01 


69.39 


-118.43 


20.08 


149.94 


4.48 


ILC, KpO 


-.439 


41.26 


63.49 


-115.40 


24.60 


150.32 


9.25 


ILC 


-0.386 


-82.66 


68.56 


66.02 


18.54 


-20.46 


-10.38 


KpO 


0.771 


-159.78 


83.95 


172.41 


-5.35 


82.67 


2.80 


Dipole 




-96.15 


48.25 










CTH (Virgo) 




-110. 


60. 










w (2,l,2) 




-105.73 


56.62 











the issues of non-Gaussianity and anisotropy. In fact there 
are three independent aspects to quadrupole anomalies: the 
low value of the angular power spectrum (C2), its possible 
non-Gaussianity, and its possible anisotropy. 

An alternative approach is that pioneered in iMagueiid 

lll995l) (the two approaches actually coincide for 1 = 1). By 
rewriting the spherical harmonics in Cartesian coordinates 
(instead of their more popular rendition in polar coordi- 
nates) one obtains a symmetric multilinear form of order I. 
Thus the i multipole may be regarded as a traceless, sym- 
metric, rank i tensor. From this one may extract 21 — 2 
(invariant) contractions, and a orthonormal frame. For ex- 
ample the quadrupole may be regarded as a 3 x 3 symmet- 
ric traceless matrix. By considering Cartesians x r = (z, x, y) 
(with constraint x 2 + y 2 + z 2 = 1) one has: 



3V6detQ 



Y. 



±1 



Y. 



±2 




y ) 



{x ± iy)z 



^-(x±iy) 2 



Thus the quadrupole may be rewritten as 

Yzm = QijX l x 3 



ST2 = S 0,2 



(9) 



(10) 



(11) 



(12) 



(15) 



if we want to separate power spectrum effects from non- 
Gaussianity issues. The variables C2, If, an d the vector sys- 
tem are now independent random variables; C2 has a xi 
distribution, if a uniform distribution in [—1,1], and the 
eigenvectors are uniformly distributed. 

The issues of isotropy and non-Gaussianity are now 
fully separated. We have checked with Monte Carlo simu- 
lations that indeed 7| and the eigenvectors are uniformly 
distributed when considering a Gaussian as seen through 
WMAP. Do the quadrupole eigenvectors point in a suspi- 
cious direction? We plot results in Table for the inverse 
noise squared coadded map, with KpO mask, for the inter- 
nal linear combination map with and without the KpO mask, 
and for the mask alone. None of the eigenvectors appears to 
correlate with th e kinematic dipole dir ection, with the direc- 
tio ns spotted by ISchwa rz&_a]| l|2004h . or with those found 
by lOliveira-Costa fc all J2004I) (which is roughly in the di- 
rection of the Virgo cluster). 

We conclude that the bispectrum analysis provided 
in this paper does not show anything abnormal for the 
quadrupole. 



leading to matrix 



Qi. 




— V 2a 22 




The 5 d.o.f. in this matrix split into 2 invariants and the 3 
d.o.f. associated with coordinate axes' orientation; for exam- 
ple 2 independent eigenvalues (the 3 eigenvalues must add 
to zero) and 3 orthonormal eigenvectors. The two invariants 
may also be the square of the matrix (the sum of the squares 
of its eigenvalues) and its determinant (the product of the 
eigenvalues). The first invariant is the power spectrum, the 
latter the bispectrum. However for a Gaussian these are not 
statistically independent, since in 3D they must satisfy the 
constraint 

13 



\detQ\ 



3^6 



(14) 



Thus, one should consider the normalised bispectrum, here 
rewritten as: 



6 DISCUSSION: SELECTION EFFECTS AND 
"ANOMALIES" 

In the bispectrum studies above we were very conservative 
about selection effects. For instance, the flat patch found in 
Section El was dismissed on the grounds that any Gaussian 
realization will reveal a flat patch somewhere if a sufficiently 
long string of is is considered (a similar argument did not 
degrade the detection in Section We stress that the flat 
patch effect may be real if we take a closer look at selection 
effects, the meaning of anomalies and the role of priors. 

"Anomaly": the very name conjures images of preju- 
dice. And yet, it's extremely difficult to make statistical 
statements without an a priori expectation. This expecta- 
tion defines the "non- anomalous" model. But if the model 
makes a large number of statistical predictions, looking into 
a sufficient number of observations will eventually turn out 
an anomaly even if the model is right. Using this argument 
we rejected the flat patch described in Section |3] However, 
it could also be that the anomaly is real, and is in fact a 
vindication of an alternative model. In this case our choice 
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of prior and subsequent classification of the anomaly as a 
fluke has dramatically prejudiced the analysis. 

An interesting historical examp le is the flat patc h found 
in the COBE-DMR four year data (lMagueiidl2000t) . Taken 
on its own it could not be dismissed as a selection effect (it 
covered most of the range to which the experiment was sen- 
sitive). If WMAP had confirmed a flat patch in the COBE 
range [i = 4 — 20) but not elsewhere, then we would have 
the strange result that more data would decrease the sig- 
nificance of a detection. This obviously means that we have 
excessively favoured a Gaussian prior. 

More interesting is what actually happened: WMAP 
failed to confirm the COBE flat patch. This showed that 
indeed the obser ved COBE flat patch was n ot a fluke: it 
was a systematic. iMaeueiio fc Medeirosl i2003l) explained in 
detail the origin of the systematic error. 

The rejection of the flat patch described in Section[3]as 
a selection effect could be a further example of the dangers 
of blind use of a Gaussian prior. Let us consider a simplified 
model in which N is the number of observed independent 
bispectrum components and p be the probability that each 
one of them is within n sigma of the mean. If we define a flat 
patch with n = 1 we have p « 0.68. The probability of a flat 
patch longer than L starting at a given fixed multipole is p L . 
In WMAP a flat patch occurs with length L — 16, starting 
at I — 32. Its probability for a Gaussian is pl = 0.002, 
seemingly ruling out Gaussianity. 

In Section [3] we argued that a selection effect was at 
play. Using our simplified model, the probability of a Gaus- 
sian generating a flat patch longer than L starting anywhere 
is p L (N — L + 1). In this case, considering that N w 100 
(number of signal dominated bispectrum components), we 
get a probability 0.18, i.e. consistency with Gaussianity. In 
Section |H] we refined this argument with Monte Carlo simu- 
lations. 

The problem with this argument is that it is reversed 
if there is a theory predicting a flat patch in the range £ = 
32 — 62, and this is used as a prior. Then Gaussianity is 
ruled out at the 99.6% confidence level. So the prior is all 
important. As the example of COBE-DMR shows the threat 
of systematics should make us beware of simply dismissing 
anomalies as flukes. 

We stress that none of these remarks apply to the re- 
sults in Section 14.11 concerning North-South asymmetries 
in the normalised inteT-£ bispectrum. We were particularly 
careful to compare like with like - letting the Gaussian sim- 
ulations choose the axis of maximal symmetry, just as we 
did with the data. The fact that the maximal asymmetry 
axis correlates well with that found for the power spectrum 
asymmetry reinforces the detection. For a Gaussian process 
power spectrum and normalised bispectrum are independent 
random variables: one would need two independent flukes 
instead of one. 



7 SUMMARY 

In this paper we have performed tests of Gaussianity and ex- 
plored possible asymmetries of the WMAP first year data. 
We calculated the inter-i?, jf, and single-^, if, normalised 
bispectrum and compared the WMAP results with those 
from simulations. We found the whole sky results to be con- 



sistent with Gaussianity, barring a 'flat patch' anomaly in 
the if for I = 32 - 62. Depending on our attitude toward 
priors this may or may not be seen as rejecting Gaussianity 
at the 99.6% confidence level. 

We also searched for hemisphere asymmetries in the bis- 
pectrum, and using the Jf found a North-South asymmetry, 
with significance level maximised at 99% for a North pole at 
(£,b) ~(57,10). This is the same North pole that maximises 
angular-power spectrum asymm etry jEriksen fc all l2004a; 
lHansen. Bandav fc G orski 20041). A similar analysis with if 
showed consistency with Gaussianity. 

We noted t hat asymmetries in the 3-point correlation 
function, C 3 (0), feriksen fc al2004al) do not diminish when 
one excludes low-^ contributions. However such asymme- 
tries are unlikely to be related to those reported in this 
paper. Our quadrupole analysis failed to find any correla- 
tion with previously repor ted anomalies dSchwarz fc al2004l : 
lOliveira-Costa fc aJl20oll 

Whether these features are real or the hallmark of fur- 
ther systematic errors remains to be seen. 
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